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Abstract 

Non-relativistic potential models are considered of the pure power V(r) = sgn(q)r q 
and logarithmic V(r) = ln(r) types. Envelope representations and kinetic potentials 
are employed to show that these potentials are actually in a single family. The log 
spectra can be obtained from the power spectra by the limit q — > taken in a smooth 
representation P n e{q) for the eigenvalues E ne (q). A simple approximation formula is 
developed which yields the first thirty eigenvalues with error < 0.04%. Extensions to 
potentials with linear combinations of terms such as -a/r + br and applications to 
spatially-symmetric few-body problems are discussed. 

PACS 12.39.Pn; 12.39.Jh; 03.65. Ge. 
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1. Introduction 

We consider first a single particle that moves in a central potential V(r) and 
obeys non-relativistic quantum mechanics. We study two cases: (a) V(r) = sgn(q)r q , 
and (b) V(r) = ln(r). We shall show that these two problems are intimately related 
and that this allows us to construct a single formula which yields accurate approx- 
imations for the corresponding discrete Schrodinger eigenvalues. The power-law and 
logarithmic potentials continue to be employed as non-relativistic models for quark 
confinement [1-7]. Hence it is important to understand that they are, in a sense, from 
the same family of potentials. A detailed analysis of the properties of non-relativistic 
potential models is given in the review article by Quigg and Rosner [8] . The approach 
to this question based on the expression {dr q /dq} q ^ = lnr is also been discussed in 
Ref. [9]. In the present paper we approach the same problem by using what we call 
the ' P -representation' for the discrete spectrum generated by power-law potentials. 
This allows us to construct a general energy eigenvalue formula which joins q = -1 
smoothly to q = 2, passing through the logarithmic case q = 0. 

The power-law family and logarithmic potential will be kept distinct until it 
becomes clear in what sense the logarithmic potential corresponds to the limit q — » 
of the power-law class. We first treat the elementary issue of scaling. We write the 
eigenvalues of the bare problems as: 

-A + sgn(q)ri — ► E nl {q) 

(1.1) 

-A + ln(r) 

where n = 1,2,3, ... counts the discrete eigenvalues in each angular-momentum sub- 
space labelled by 1 = 0,1,2, The eigenvalues so labelled have degeneracy precisely 

2£ + 1. Scaling arguments [10,3] can be used to show that more general Hamiltonians 
have the corresponding eigenvalues given in terms of {E ni (q), E% e } by 

-wA + v sgn(q)r* — to (£) & E nt {q) 

(1.2) 

-uA + vln(r) v E^ e - Ivln(^) , 

where v > is a coupling parameter and u > might, for example, be w = h 2 /(2m). 
As a consequence of these scaling rules we shall only need to consider the special case 
lu = v = 1 for this family of potentials. 

In Section (2) we review the concepts of envelope representations and kinetic 
potentials [10, 11]. These geometrical methods allow us to construct an exact semi- 
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classical representation for the discrete spectra of Schrodinger operators in which 
the eigenvalues are recovered by a minimization over a single real kinetic-energy 
variable. In Section (3) we show, for the power-law family of potentials, that this 
leads, in particular, to the smooth representation P n e{q) for the eigenvalues E(q) (with 
v = 1 ) which, in turn, allows us to show that the log spectrum can be derived as a 
limit from the family of power spectra. As a motivation for this, we plot the power 
eigenvalues E nl (q) as functions of q in Fig.(l). These graphs are to be compared 
with the corresponding smooth P -representation shown in Fig. (2). Once the P- 
representation is established, we go on to consider the log-power theorem P n e{0) = P^ 
and, in Section (4), to develop an approximation formula for P n e(q)- In Section (5) 
we consider potentials which are linear combinations of powers, such as the Coulomb 
plus linear potential -a/r + br; in Section (6) we use these results to study few-body 
problems for which spatially-symmetric states are accessible. 

2. Envelope representations and kinetic potentials 

We distinguish a potential V(r) = vf(r) from its shape /(r), where the positive 
parameter v is often called the 'coupling constant'. The idea behind envelope repre- 
sentations is suggested by the question: if one potential (shape) f(r) can be written 
as a smooth transformation /(r) = g{h(r)) of another h(r), what spectral relationship 
might this induce? We consider potential shapes that support at least one discrete 
eigenvalue for sufficiently large values of the coupling v and suppose for definiteness 
that the lowest eigenvalue of -A + vh(r) is given by H(v) and that of -A + vf(r) 
by F(v). If the transformation function g(h) is smooth, then each tangent to g is 
an affine transformation of the 'envelope basis' h of the form /*(r) = a(t)h(r) + /3(t), 
where r = t is the point of contact. The coefficients a(t) and j3(t) are obtained by 
demanding that the 'tangential potential' f^\r) and its derivative agree with f(r) 
at the point of contact r = t. Thus we have 

a(t) = f'(t)/h'(t), = f(t) - a(t)h(t). (2.1) 

The geometrical configuration is illustrated in Fig. (3) in which f(r) = -1/r + r and 
the envelope basis for the upper family is the linear potential h{r) = r and, for 
the lower family, the Coulomb potential h(r) = -l/r. The spectral function for the 
tangential potential f {t \r) = a(t)h(r) + (3{t) is given by H^(v) = H(va(t)) + v(3{t). If 
the transformation g(h) has definite convexity, say g"(h) > 0, then each tangential 
potential lies beneath f(r) and, as a consequence of the comparison theorem, we 
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know that each corresponding tangential spectral function and the envelope 

of this set, lie beneath F(v). Similarly, in the case where g is concave ( g"(h) < 0), we 
obtain upper bounds to F(v). These purely geometrical arguments, depending on the 
spectral comparison theorem, extend easily to the excited states. The spectral curves 
corresponding to the envelope representations for the potential in Fig. (3) are shown 
in Fig. (4) for the excited state (n,£) = (2,4). For comparison the exact curve E = F(v) 
is also shown; this curve will be close to the Coulomb envelope for large v and to the 
linear envelope for small v. Of course, the envelopes of which we speak still have to 
be determined explicitly. 

Whilst extensions of this idea to completely new problems, such as simultaneous 
transformations of each of a number of potential terms [12], are best formulated 
initially with the basic argument outlined above, the question arises, in the 1-term 
case, as to whether there is a simple way of determining the envelopes of the families 
of upper and lower spectral functions. One solution of this problem is by the use 
of 'kinetic potentials' which were introduced [10, 11] precisely for this purpose. The 
idea is as follows. To each spectral function F ne (v) there is a corresponding 'kinetic 
potential' (= minimum mean iso-kinetic potential) f n e{s)- The relationship between 
F n £ and f n e is invertible and is essentially that of a Legendre transformation [14]: we 
can prove in general that F is concave, f(s) is convex, and 

*»/"(*) = -1. (2.2) 

The explicit transformation formulas are as follows 

{f ni (s) = F^ e (v), s = F ne (v)-vF^ e (v)} <- {F ni (v)/v = f ni (s)-sfM, l/v = -f' nl {s)}. 

(2.3) 

An a priori definition of the ground-state kinetic potential fw(s) = f(s) is given by 

f(s)= inf (2-4) 
(V,V)=i 

(V>,-Ai/>)=s 

where V c L 2 (IR 3 ) is the domain of the Hamiltonian. The definition for the excited 
states is a little more complicated [11] and, in view of (2.3), will not be needed in 
what follows. 

What is important is that the spectral functions, either exact or approximate, 
are recovered from the corresponding kinetic potentials by a minimization over the 
kinetic-energy variable s. The total minimization required by min-max has been 
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divided into two steps: the first is constrained by (tp, -Aip) = s and the second is a 
minimization over s. We have in all cases: 



Fnt(v) = mm {s + vf nl (s)} . 



(2.5) 



Another form of this expression is possible for the kinetic potential is monotone and 
allows us to change variables (s — > r) by /(r) = f n e(s). Thus we have 

F nl {v) = mm {K«\r) + vf(r)} , K<£ = /"/ o /. (2.6) 

The two corresponding expressions of the envelope approximation then become 

U{s)^f{K e {s)), or KflnKM. (2.7) 

The second form (2.6) of the expression for F ni (v) isolates the potential shape / itself 
and leads to an inversion sequence [15] which reconstructs the potential from a single 
given spectral function; but this is another story. 

It is useful here to provide the formulas for the kinetic potentials corresponding 
to pure power-law potentials. Since, by elementary scaling arguments, we have 



-A + v sgn(q)r q 



F(v) = F(l)v~m. 



(2.8) 



where = E(q) is the bottom of the spectrum of -A + sgn(q)r q , we immediately 
find from (2.3) that the kinetic potentials for these potentials are given by 



fni(s) 






qE ni (q) 




2 + q 



s 2 . 



Meanwhile the corresponding if -functions all have the same simple form 

*#(r)=(^y, 

where the P numbers are given by 



Pne(q) - \E ni (q)\ 2 -^ 



2 + q 



q^O. 



(2.9) 



(2.10) 



(2.H) 



Consequently the power-law kinetic potentials (2.9) may be expressed in the simple 
form 

'Pndq)Y 



fni(s) = sgn(q) 



S2 



q^O. 



(2.12) 



Meanwhile, as we shall see, for the log potential f(r) = log(r) we have f„e(s) = 
log(P n e(0)/si). These formulas are particularly useful for the analysis of sums of powers 
discussed in Section (5) below. 
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3. The functions E(q) and P(q) for pure powers and the log potential 

By using the P -representation for the spectrum we have from (2.6) and (2.10) 

^ = ^{^ + y(r) }' (3 ' 1} 
Independently of the theory discussed in Section (2), we see that this representation 
is well defined provided we can always construct the P which corresponds to a given 
E. For smooth potentials, the minimum exists and is unique if r 3 V'(r) is monotone, a 
condition which is certainly satisfied if V(r) is a pure power or log. For these poten- 
tials, Eq.(3.1) is therefore an exact semi-classical representation for the eigenvalues: 
the kinetic energy term scales like L~ 2 , as it should. It is perhaps interesting also 
to note that positive factors such as u> and v can be re-introduced in front of the 
kinetic- and potential-energy terms, without having to revise the P numbers. For the 
specific pure-power and log problems at hand, we have respectively: 



E ne (q)= mm \^4^+sgn(qy 



r>0 



= sgn(q) (; 



1 , ,\ f 2P nM\« + 



(3.2) 



(3.3) 



,2 

and 

= ln((2e)*P£). 

Some of the E, and the corresponding P, are already known exactly, from elemen- 
tary quantum mechanics. Thus from the known eigenvalues for the Hydrogen atom 
Ent(—1) = -[2(n + l)]~ 2 and the harmonic oscillator E ne (2) = 4n + 2£- 1 we immediately 
obtain the outer P numbers of our range of interest -1 < q < 2 -. 

P nt (-l)=n + £, P nt (2) = 2n + £-^. (3.4) 

Our goal is now to interpolate P n e(q) at interior points of the range. First we use 
kinetic potentials to establish two theoretical results. It is straightforward to show 
that each power potential with power q x is a convex transformation of all other power 
potentials with powers q 2 < qi- Meanwhile, the log potential is a convex function of 
each negative power, and, at the same time, it is a concave transformation of each 
positive power. Two results which follow from this analysis are (i) [13] the functions 
P n e(q) are monotone increasing in q, and (ii) [16] the log-power theorem, which says 
that lim^o Pra(q) = PndO) = 
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As is clear from Fig.(l), the spectra E n i(q) of the power potentials converge on 
the values ±1 as the power q approaches zero respectively from positive or negative 
values. Meanwhile, the set of points for the log spectrum would appear to be quite 
unrelated, and would have no natural place on this graph. In the P -representation the 
picture is quite different, as is shown in Fig. (2). The monotone P(q) curves are distinct 
at q = 0, and we know that they are exactly equal there to the values corresponding 
to the log potential. This smooth expression of the entire spectral family suggests that 
we can interpolate the P curves with a simple polynomial, a task which is taken up 
in the next Section. We made an earlier attempt at an interpolation [13], which was 
much less accurate because it could not benefit from the (yet undiscovered) log-power 
theorem. 

4. Interpolation formula for P{q) 

As a model for the P curves we use a cubic polynomial which allows us to fit values 
at the four points q = {-1,0,1,2}; the exact formulas (3.4) are used for q = {-1,2}. 
We make the expansion about q = so as to favour the q values between log and 
linear that are used in potential models. Thus we define 



P(q) =a + bq + cq 2 + dq 3 



(4.1) 



in which, for simplicity, we have omitted the quantum-number subscripts nt. By 
inversion we have 



"a" 




b 




c 













1 





1 

3 

2 





_ 1 

6 



1 

6 





-P(-1)1 




P(0) 




P(l) 




L P(2) J 



(4.2). 



For the linear potential q = 1 the exact S-state eigenvalues may be expressed in terms 
of the zeros of the Airy function [17]. We take these known values and complement 
them with others computed numerically for the log and linear potentials. By inverting 
Eq.(3.2) and Eq.(3.3) we have thus computed the coefficients P(0) for log and P(l) 
for linear. The P-data for the first 30 eigenvalues are exhibited in Table(l), along with 
the approximations we get for the energies E n i(\). For comparison, we have tabulated 
also the percentage errors in these eigenvalue results. For the first 30 eigenvalues, the 
percentage errors are all positive and less than 0.04% at q = \; as q approaches 
q = or q = 1, these errors decrease dramatically. Similar accuracy is obtained for 
the whole range -1 < q < 2 of the interpolation. 
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5. Linear combinations 

The kinetic-potential formalism turns out to be helpful for another class of poten- 
tial composition, namely linear combinations. Consider the bottom of the spectrum 
of the problem with hamiltonian 

H = -A + af (1) (r)+bf (2 \r). (5.1) 

Suppose that V is the exact normalized ground-state wave function corresponding to 
eigenvalue E and that < u < 1, then we have exactly 

H = u> (-A + ^/«(r)) + (1 - w) (-A + y^-Z^m) 

and consequently for E = (ip,Hip) we find the lower bound 

E> ajF ^(-) + (l-w)F ( - 2 \— — ). (5.2) 
u> 1 — u> 

By maximizing the right-hand side of (5.2) with respect to u> we obtain a best lower 
bound to E expressed in terms of the 'component' eigenvalues. It turns out [11] 
that this lower bound (after optimization over u) is given exactly by applying the 
kinetic-potential rule 

f(s)>ap\s) + bfW(s). (5.3) 

This lower bound has been extended to arbitrary linear combinations and to con- 
tinuous mixtures [18]; from the argument sketched above it is clear that the same 
bound is valid for the bottom of each angular-momentum subspace. What is perhaps 
surprising is that for the other excited states it remains a 'good' approximation (error 
a few %). It is in effect an approximate generalization of the classical operator sum 
theorem of Weyl [19- 21]. The recipe we obtain for the eigenvalues of the Hamiltonian 
H = -A + J2 q a(q)sgn(q)r q for a(q) > becomes, from (2.5) and (2.12), with the change 
of variable s = l/r 2 , as follows 

l + ]Ta(g)sgn(q) {P nl (q)r) q 
i 

This expression is an energy lower bound for the bottom of each angular-momentum 
subspace n = 1. Thus kinetic potentials are almost additive; for the bottom of each 
angular-momentum subspace they are sub-additive; whenever the sum has only one 
term, the result is exact. 

As an example we consider the Coulomb plus linear potential f(r) = -a/r + br and 

find 

r 1 n i 

(5.5) 



E nt w min 

r>0 



(5.4) 



E nt w min 

r>0 



1 a L 

— V bfir 

r z vr 
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where we have a lower bound for Coulomb envelopes v = fj, = P n e(-1) = n+£, we have an 
upper bound with linear envelopes v = ^ = P n e(l), and we obtain an approximation by 
using the sum approximation v = P n e(-1), ^ = P n tiX), which latter is a lower bound 
whenever n = 1. It is interesting that an upper bound to the ground-state energy E w 
derived from a scale optimized Gaussian wave function is also provided by the same 
formula (5.5) if we assign 

37r\ * / 6 \ * 

— = 1.085401, n = - = 1.381977. (5.6) 



The Gaussian upper bound will be important for the TV-body problem discussed in 
the next section. 

Some progress can even be made towards a formula for these approximations [26]. 
By scaling arguments we have for every discrete eigenvalue £(w,a,b) of 

2 i. 2 

H = -wA - - + br, -» £(w, a, b) = —£ (1, 1, A) , A = (5.7) 

By applying (5.5) with a = 1 and b = A we obtain an expression for E(X) whose 
inverse can be written explicitly as 

= 2(^)3 -^ 2 [(1 + 3^)^-1] 
M [(l + 3^ 2 S)i -l] 3 

where the values of v and n for the various bounds and approximations are given 
above, and E > -l/4v 2 . This is a global formula valid for all positive a and b. 

6. TV-body problems 

We now discuss an application to TV-body problems for which a spatially- 
symmetric state is accessible, such as a simple non-relativistic model for a 3-quark 
system representing a nucleon. Some general features and results for specific many- 
body systems may be found in Refs.[22, 23]. The Hamiltonian, with center-of-mass 
removed, for a system of N identical particles each of mass m interacting via central 
pair potentials may be written 

1=1 \i=l / 3>i=l v 7 

where V a and a are respectively the depth and range parameters of the potential with 
shape /. By algebraic rearrangement (6.1) may be rewritten in the more symmetrical 
form 



j>i=l 
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We now define new coordinates by p = BR, where p = [p»] and R = [r,] are column 
vectors of the new and old coordinates, respectively, and B is a real constant N x N 
matrix. For convenience we require all the rows of B to be unit vectors, we let the 
elements of the first row all be equal to -^=, so that p\ is proportional to the centre- 
of-mass coordinate; we also require that the remaining N—l rows of B be orthogonal 
to the first row, so that they define a set of N—l relative coordinates. One more row 
is also fixed so that we have at least one pair distance at our disposal, namely 

P2 = ^r- (6 - 3) 

For spatially-symmetric states, we have found that Jacobi relative coordinates, for 
which B is orthogonal, are the most useful. Thus, corresponding to the transformation 
p = BR of the coordinates, it follows that the column vector P of the associated 
momenta transforms to the new momenta n = [7^] by the relation n = (B T )~ 1 P = P. 
If * is any translation-invariant wave function for the TV-body system composed 
of identical bosons, then we can write [27, 28] the following mean energy relation 
between the iV-body and 2-body systems: 

(*,W*) = (*,H*) ) (6.4) 

where the 'reduced' two-particle Hamiltonian IH is given by 

Further simplifications can be achieved if we work with dimensionless quantities. 
We suppose that the translation-invariant TV-body energy is £ and we define the 
dimensionless energy and coupling parameters E and v by the expressions 

m£a 2 NmV a 2 

E= tt, v = o — . 6.6 

{N-l)h 2 2H 2 V ' 

It is then natural to define a dimensionless version of the reduced 2-body Hamiltonian 
IH and the relative coordinate p 2 by the relations 

g= (S = " A + " /(r)l r = V2p 2 /a = r 1 -r 2 , r=||r||. (6.7) 

We note that the Hamiltonian H depends on N only through the dimensionless 
coupling parameter v. By the Rayleigh-Ritz (min-max) principle [29-31], we have the 
following characterization of the TV-body ground-state energy parameter E in terms 
of H : 

E = u i n^ ; ^=F N (v), (6.8) 
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where * is a translation-invariant function of the N - 1 relative coordinates (and 
spin variables, if any) which is symmetric or antisymmetric under the permutation 
of the N individual-particle indices. The iV-body energy £ is recovered from E by 
inverting (6.6). Thus we have explicitly: 

S = {N ~? 2 F N (>™£) . (6.9) 
ma 2 \ 2h 2 J K ' 

The energy bounds we use are summarized in terms of the F functions 

F 2 (v) <E = F N (v)< Faoty) < F G (v). (6.10) 

The history of the equivalent 2-body method for TV-particle systems has been de- 
scribed in earlier papers [24-27] and in the references therein. The main result is a 
general energy lower bound which, with orthogonal Jacobi relative coordinates, is 
given by 

F 2 (v) <E = F N (v), (6.11) 

where F 2 (v) is the lowest eigenvalue of the 1-particle ('reduced' 2-particle) Hamilto- 
nian H = -A + vf(r). With equal simplicity, our upper bound F G (N), provided (for 
spatially-symmetric states) by a Gaussian trial function may also be expressed in 
terms of the 1-body operator H. This is a consequence of the following argument. If 
and only if the spatially-symmetric translation-invariant function * is Gaussian [32, 
33], it may be factored in the form 

*(/>2, Pa, ■ ■ ■ , Pn) = tp(p2)v(P3, Pi,---, Pn)- (6.12) 

But the equivalence (6.4) then implies, in this case, that E < (ip,H?p)\\'ip\\-' 2 . This 
explains why the inequalities (6.10) collapse together in the case of the harmonic os- 
cillator to the common value E = 3vi for which the exact 1-body lowest eigenfunction 
of H is also Gaussian. In this argument we assume for the upper bound that <H> 
has been optimized with respect to the scale of the wave function. The symmetry 
of these Gaussian functions is demonstrated most clearly by the following algebraic 
identity: 

N N 

£ (iW^A^p^ (613) 

j>i=l k=2 

As an example [34] we consider a model symmetric TV-quark system and we take 
some physical numbers from two early papers by Kang and Schnitzer [35], and Gromes 
and Stamatescu [36], and allow 2 < N < 5 for this illustration, although N > 3 no 
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longer corresponds to a physical system. The units are % = c = 1 and we assume equal 
masses m = 0.3GeV and Coulomb-like coupling a = 0.35. We keep a fixed and plot 
the energy obtained from the formula (5.5) with v = 1 and \i = 2|i?i (l)/3| 3 / 2 = 1.376083 
for a lower bound, and the same formula with the Gaussian values (5.6) for an upper 
bound. The energy curves as a function of the linear coupling b are shown for each 
N in Figs. (5) and (6), which only differ in the range of b. One can only expect to 
get modest benefit from such a simple model. However, definite energy bounds could 
provide a small island of security in an otherwise highly complicated many-body 
environment. 

7. Conclusion 

In spite of present-day computing convenience, it is still very useful to have an 
approximate eigenvalue formula, particularly a simple and accurate one. The possi- 
bility of a simple formula is a consequence of the existence of the smooth monotone 
P -representation for the eigenvalues. Since the early work on potential models there 
has been a realization that, from a practical point of view, a log and a power poten- 
tial, with < q < 1, could serve in the model almost equally well, especially if q is 
small. The log-power theorem provides a theoretical basis for these concrete spectral 
observations. Once we have the smooth P representation for the power-law and log 
eigenvalues we are able to provide spectral information concerning linear combina- 
tions of such terms and also for ./V-body problems such as non-relativistic models for 
systems of few quarks. 
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Table 1 The 'input' values of P n e(0) and P n e(l); and the approximations E^ e {\) for 
E n t{\) obtained via the cubic P formula (3.1), with the percentage errors. 



n 


n 
I 




PntW 


rpA ( 1 \ 
E n t\2> 


771 / 1 \ 


0/ 
/o 


1 


U 


1 Ol 0£?7 

l.zlob / 


l.o (buo 


l.OOO/O 


l.ooooy 


U.Uly 


Z 


U 


Z. 1 zUbO 


Q 1Q1Q1 

0. lolol 


O c; c; 1 /io 
Z.0014Z 


z.OOUbO 


U.UoU 


Q 

o 


U 


4.ZOO00 


4.yyzoo 


O.U0ZZ4 


q nc;i 1 q 
o.UOllo 


U.UoO 


A 

4 


U 


0. /4ybz 


b.oU014 


Q /I £Q /1 1 

0.40o41 


Q /I £01 Q 

o.40zlo 


n HQ7 

U.Uo / 


r 




U 


/ .zb /Uo 


o.blozo 


Q 70/ICO 

0. f y4oZ 


0. 1 yoob 


U.Uoy 


1 


1 


Z.zlo4o 


O Q71 OO 

z.o / lyz 


onn7Q 
Z.oUU / 


Z.oUUOU 


U.U1U 


z 


1 


O.boOOO 


/lie; c;n 1 
4. loom 


O ON Of? 

Z.o04ob 


Z.o04o4 


U.Ulo 


Q 

o 


1 


0.1///4 


O.yOoUU 


o.Zoboy 


O.Z0O00 


n no £ 
U.UoO 


A 

4 


1 


o.o / yob 


/ . /0 lui 


Q RA CQ K 

o.b4oo0 


Q /?/1 70(1 

o.b4 / oy 


U.Uzb 





1 


o.lobU / 


0. c^K/ino 
y.0b4Uo 


o.ybooZ 


o.ybzbo 


n non 
U.Uzy 


1 




o oi i /in 
o.Z114y 


0.0 / Ulo 


Z.bOf fo 


Z.bO /0b 


n nn7 
U.UU I 


Z 




4.00O0U 


£ 1 /1 1 qc: 
0. 141o0 


o.lZU 1 1 


0. IzUoo 


n m /i 
U.U14 


Q 
O 


Z 


0.14b / Z 


b.yzyl 1 


o.OUoUy 


o.0Uz40 


U.Ulo 


4 


Z 


/ .boboy 


Q 70 ^1 ^ 

0. 1 ZOlO 


0.00000 


o oqoc; /I 
O.OOZ04 


n noi 
U.Uzl 





z 


n 1 qq 1 
y.looly 


lU.OZOyb 


/I 1 OK7Q 

4.1Zb /o 


/| 1 ticoi 
4. IZOol 


n no/i 
U.UZ4 


1 


Q 
O 


a 01 n/i a 
4.Z1U44 


/I Q(?OOQ 

4.obyZo 


Z.y04bl 


ncc A A £ 
Z.90440 


U.UUO 


Z 


Q 
O 


O.bOo I'd 


b.loZyo 


0.0O / yo 


Q Q£7£n 

0.0O /0y 


n m 
U.Ulz 








7 1 Ofififi 

( .izbob 


7 Ol QD/l 

( .yiou4 


Q 7DQ07 
0. ( UOZ 1 


0. IvZ l\J 


U.Ulo 


4 


3 


8.60714 


9.70236 


4.00810 


4.00737 


0.018 


5 


3 


10.09555 


11.49748 


4.28283 


4.28196 


0.020 


1 


4 


5.20980 


5.36863 


3.21247 


3.21233 


0.004 


2 


4 


6.65235 


7.12732 


3.57310 


3.57275 


0.010 


3 


4 


8.11305 


8.90148 


3.88950 


3.88898 


0.013 


4 


4 


9.58587 


10.68521 


4.17335 


4.17268 


0.016 


5 


4 


11.06163 


12.47532 


4.43164 


4.43131 


0.008 



Smooth spectral transition from Coulomb to oscillator 



page 16 




Figure 1. 

The first 30 eigenvalues E nt {q), 1 < n < 5, < ^ < 5, corresponding to the power 
potential V(r) = sgn(q)r q . For q > 0, the eigenvalues increase with q from 1 to 
E nt {2) = 4n+2f+l; for q < 0, they decrease (as q increases) from E ne (-1) = -{2(n+l)}- 2 . 
to -1. Both sets of curves increase with n and I. 
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Figure 2. 

In the P -representation, the same set of 30 eigenvalues shown in Fig.(l) now lie on 
monotone smooth curves. The log-power theorem states that the P values for the 
log potential are precisely P n e(0)- As q increases from -1 to 2, the degeneracy of 
the Coulomb problem P n e(-1) = n + £ evolves into the degeneracy of the harmonic 
oscillator P ne (2) = 2n + (.-\. 
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f(r) = -a/r + br 



f(r) 




b = 1 



Figure 3. 

The linear plus Coulomb potential shape f(r) = -1/r + r represented as the envelope 
curve of two distinct families of potentials of the form ah(r) + (3. In the upper family 
h{r) = r is a linear potential; in the lower family h(r) = -1/r is a Coulomb potential. 
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V(r) = v(-a/r + br) 




1 2 3 4 5 



Figure 4. 

The spectral approximation corresponding to Fig. (3). Each 'tangential' potential 
/W(r) = ah(r) + (3 generates a corresponding tangential energy curve F^\v) = H(av) + 
(3v. The envelopes of these spectral families generate upper and lower bounds to the 
exact curve E = F(v), shown here for the case (n,£) = (2,4). 
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V(r) = -0.35/r + br 
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Figure 5. 

In a naive model with linear plus Coulomb pair potentials, the graphs show upper and 
lower bounds to the energies of the lowest spatially many-body symmetric iV-quark 
states. These states are only physically accessible for N < 3, but the bounds are valid 
for all N > 2. In the case of TV = 2 there is only one curve because the solution is 
'exact'. 
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V(r) = -0.35/r + br 
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Figure 6. 

The many-body graphs are the same as in Fig. (5) but with a smaller scale showing 
the intersection of the curves. Graphs such as this may help with the elementary 
consistency question: can we devise a potential that yields the known 2- and 3-quark 

masses? 



